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We propose a general mechanism of random-field-induced order (RFIO), in whicii long-range order is induced 
by a random field that breaks the continuous symmetry of the model. We particularly focus on the case of the 
classical ferromagnetic XY model on a 2D lattice, in a uniaxial random field. We prove rigorously that the 
system has spontaneous magnetization at temperature T = 0, and we present strong evidence that this is also 
the case for small T > 0. We discuss generalizations of this mechanism to various classical and quantum 
systems. In addition, we propose possible realizations of the RFIO mechanism, using ultracold atoms in an 
optical lattice. Our results shed new light on controversies in existing literature, and open a way to realize RFIO 
with ultracold atomic systems. 
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I. INTRODUCTION 



A. Disordered ultracold quantum gases 



Studies of disordered systems constitute a new, rapidly de- 
veloping, branch of the physics of ultracold gases. In con- 
densed matter physics (CM), the role of quenched (i.e. in- 
dependent of time) disorder cannot be overestimated: it is 
present in nearly all CM systems, and leads to numerous phe- 
nomena that dramatically change both qualitative and quanti- 
tave behaviours of these systems. This leads, for instance, to 
novel thermodynamical and quantum phases'-^, and to strong 
phenomena, such as Anderson localization^'*-^-^. In general, 
disorder can hardly be controlled in CM systems. In con- 
trast, it has been proposed recently, that quenched disorder 
(or pseudo-disorder) can be introduced in a controlled way in 
ultracold atomic systems, using optical potentials generated 
by speckle radiations^^, impurity atoms serving as random 
scatterers"^, or quasi-cristalline lattices". This opens fantastic 
possibilities to investigate the effect of disorder in controlled 
systems (for a review in the context of cold gases, see Ref.[l2h. 
Recently, several groups have initiated the experimental study 
of disorder with Bose-Einstein condensates (BEC)'^''*'^, and 
strongly correlated Bose gases'^. In the center of interest of 
these works is one of the most fundamental issues of disor- 
dered systems that concerns the interplay between Anderson 
localization and interactions in many body Fermi or Bose sys- 
tems at low temperatures. In non-interacting atomic sytems, 
localization is feasible experimentalljiii, but even weak inter- 
actions can drastically change the scenario. While weak re- 
pulsive interactions tend to delocalize, strong ones in confined 
geometries lead to Wigner-Mott-Uke locaUzationi^. Both ex- 
periments and theory indicate that in gaseous systems with 
large interactions, stronger localization effects occur in the 
excitations of a BEGi^ii2i22iSi,, rather than on the BEC wave- 
function itself. In the limit of weak interactions, a Bose gase 



enters a Lifshits glass phase, in which several BECs in vari- 
ous localized single atom orbitals from the low energy tail of 
the spectrum coexist^^ (for 'traces' of the Lifshits glass in the 
meanfield theory, see Ref. .15.) . Finally, note that disorder in 
Fermi gases, or in Femi-Bose atomic mixtures, should allow 
one to realize various fermionic disorderded phases, such as 
a Fermi glass, a Mott-Wigner glass, 'dirty' superconductors, 
etc. (Ref.ll2l). or even quantum spin glasses 
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B. Large effects by small disorder 



One of the most appealing effects of disorder is that even 
extremely small randomness can have dramatic consequences. 
The paradigm example in classical physics is the Ising model 
for which an arbitrarily small random magnetic field destroys 
magnetization even at temperature T=0 in two dimensions, 2D 
(Refs. 124 25), but not in D > 2 (Ref. IH). This result has 
been generalized to systems with continuous symmetry in ran- 
dom fields distributed in accordance with this symmetry^i2^. 
For instance, the Heisenberg model in a S'0(3)-symmetrically 
distributed field does not magnetize up to 4D. 

In quantum physics, the paradigm example of large ef- 
fects induced by small disorder is provided by the above- 
mentioned Anderson localization which occurs in ID and 2D 
in arbitrarily small random potentials^. In this paper, we pro- 
pose an even more intriguing opposite effect, where disor- 
der counter-intuitively favors ordering: a general mechanism 
of random- field- induced order (RFIO) by which certain spin 
models magnetize at a higher temperature in the presence of 
arbitrarily small disorder than in its absence, provided that the 
disorder breaks the continuous symmetry of the system. 



C. Main results and plan of the paper 

As is well known, as a consequence of the Mermin-Wagner- 
Hohenberg theorem^^, spin or field theoretic systems with 
continuous symmetry in dimensions less or equal to 2D can- 
not exhibit long range order The mechanism that we pro- 
pose here breaks the continuous symmetry, and in this sense 
acts against the Mermin-Wagner-Hohenberg no-go rule in 2D. 
In particular, we prove rigorously that the classical XY spin 
model on a 2D lattice in a uniaxial random field magnetizes 
spontaneously at T = in the direction perpendicular to the 
magnetic field axis, and provide strong evidence that this is 
also the case at small positive temperatures. We discuss gen- 
eralizations of this mechanism to classical and quantum XY 
and Heisenberg models in 2D and 3D. In 3D, the considered 
systems do exhibit long range order at finite T > 0, but in 
this case the critical temperature decreases with the 'size' of 
the symmetry group: it is the largest for the Ising model (the 
discrete group Z2), higher for the XY model [the continu- 
ous group U{1)], and the highest for the Heisenberg model 
[the continuous group SU{2), or 5*0(3)]. In this case we ex- 
pect that our mechanism will lead to an increase of the crit- 
ical temperature for the XY and Heisenberg models, and to 
an increase of the order parameter value at a fixed temperature 
for the disordered system in comparison to the non-disordered 
one. Finally, we propose three possible and experimentally 
feasible realizations of the RFIO phenomenon using ultracold 
atoms in optical lattices. 

The paper is organized as follows. In section|lIl we present 
the results concerning the RFIO in the classical XY model on 
a 2D lattice. First, we rigorously prove that the system magne- 
tizes in the direction perpendicular to the direction of the ran- 
dom magnetic field at T = 0, and then, we present arguments 
that the magnetization persists in small T > case, as well 
as the results of numerical classical Monte Carlo simulations. 
Section |III| is fully devoted to the discussion of the general- 
izations of the RFIO mechanism to several other classical and 
quantum spin systems. In section [TVl we discuss several ex- 
perimentally feasible realizations of RFIO in ultracold atomic 
systems. Finally, we summarize our results in section IVl 



II. RFIO IN CLASSICAL XY MODEL 



A. The system under study 



Consider a classical spin system on the 2D square lattice 
Z^, in a random magnetic field, h (see Fig.[Tl). The spin vari- 
able, (7i ~ {cos 9i, sin 6i), at a site i € Z^ is a unit vector 
in the xy plane. The Hamiltonian (in units of the exchange 
energy J) is given by 



H/J = - ^ a.i ■ aj - ey^hj ■ aj 
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Here the first term is the standard nearest-neighbor interac- 
tion of the XY-model, and the second term represents a small 



Figure 1 : (color online) XY model on a 2D square lattice in a random 
magnetic field. The magnetic field is oriented along the y axis, hi = 
r\ie.y, where r\i is a real random number. Right boundary conditions 
are assumed on the outer square, possibly placed at infinity (see text). 



random field perturbation. The h^'s are assumed to be inde- 
pendent, identically distributed random, 2D vectors. 

For e = 0, the model has no spontaneous magnetization, m, 
at any positive T. This was first pointed out in Ref. (281, and 
later developed into a class of results known as the Mermin- 
Wagner-Hohenberg theorem^^ for various classical, as well as 
quantum two-dimensional spin systems with continuous sym- 
metry. In higher dimensions the system does magnetize at low 
temperatures. This follows from the spin wave analysis^'', and 
has been given a rigorous proof in Ref. 30. The impact of 
a random field term on the behavior of the model was first 
addressed in Refs. l24ll25l where it was argued that if the dis- 
tribution of the random variables hj is invariant under rota- 
tions, there is no spontaneous magnetization at any positive T 
in any dimension Z? < 4. A rigorous proof of this statement 
was given in Ref. 25. Both works use crucially the rotational 
invariance of the distribution of the random field variables. 

Here we consider the case where h^ is directed along the y- 
axis: h; = r\ie,y, where Gj, is the unit vector in the y direction, 
and r\i is a random real number Such a random field obviously 
breaks the continuous symmetry of the interaction and a ques- 
tion arises whether the model still has no spontaneous mag- 
netization in two dimensions. This question has been given 
contradictory answers in Refs. ! 3lll32t while Ref. [3ll predicts 
that a small random field in the jy-direction does not change 
the behavior of the model, Ref. ^ argues that it leads to the 
presence of spontaneous magnetization, m, in the direction 
perpendicular to the random field axis in low (but not arbitrar- 
ily low) temperatures. Both works use renormalization group 
analysis, with Ref. |32| starting from a version of the Imry-Ma 
scaling argument to prove that the model magnetizes at zero 
temperature. 

The same model was subsequently studied by Feldman'^-, 
using ideas similar to the argument given in the present paper 
As we argue below, however, his argument contains an essen- 
tial gap, which is filled in the present work. We first present a 
complete proof that the system indeed magnetizes at T = 0, 



and argue that the ground state magnetization is stable under 
inclusion of small thermal fluctuations. For this, we use a ver- 
sion of the Peierls contour argument^"*, eliminating first the 
possibility that Bloch walls or vortex configurations destroy 
the transition. 



B. RFIO at T = 

Let us start by a rigorous analysis of the ground state. Con- 
sider the system in a square A with the 'right' boundary con- 
ditions, ai = (1, 0), for the sites i on the outer boundary of A 
(see Fig. [T]i. The energy of any spin configuration decreases 
if we replace the x components of the spins by their absolute 
values and leave the y components unchanged. It follows that 
in the ground state, x components of all the spins are nonneg- 
ative. As the size of the system increases, we expect the x 
component of the ground state spins to decrease, since they 
feel less influence of the boundary conditions and the ground 
state value of each spin will converge. We thus obtain a well- 
defined infinite-volume ground state with the 'right' boundary 
conditions at infinity. 

We emphasize that the above convergence statement is non- 
trivial and requires a proof. Physically it is, however, quite 
natural. A similar statement has been rigorously proven for 
ground states of the Random Field Ising Model using Fortuin- 
Kasteleyn-Ginibre monotonicity techniques^^*^. 



C. Infinite volume limit 

A priori this infinite-volume ground state could coincide 
with the ground state of the Random Field Ising Model, in 
which all spins have zero x component. The following argu- 
ment shows that this is not the case. Suppose that the spin ai 
at a given site i is aligned along the y-axis, i.e. cos di ~ 0. 
Since the derivative of the energy function with respect to di 
vanishes at the minimum, we obtain 






0. 



(2) 



Since cosdi = 0, this implies J2i-\i-i\=i cos9j = 0. Because 
in the 'right' ground state all spins lie in the (closed) right 
half-plane x > 0, all terms in the above expression are non- 
negative and hence have to vanish. This means that at all the 
nearest neighbors j of the site i, the ground state spins are di- 
rected along the j;-axis as well. Repeating this argument, we 
conclude that the same holds for all spins of the infinite lat- 
tice, i.e. the ground state is the (unique) Random Field Ising 
Model ground state. This, however, leads to a contradiction, 
since assuming this, one can construct a field configuration, 
occurring with a positive probability, which forces the ground 
state spins to have nonzero x components. To achieve this 
we put strong positive {rji > 0) fields on the boundary of a 
square and strong negative fields on the boundary of a con- 
centric smaller square. If the fields are very weak in the area 
between the two boundaries, the spins will form a Bloch wall. 



rotating gradually from 6 = 7r/2 to 9 = — 7r/2. Since such 
a local field configuration occurs with a positive probability, 
the ground state cannot have zero x components everywhere, 
contrary to our assumption. 

We would like to emphasize the logical structure of the 
above argument, which proceeds indirectly assuming that the 
ground state spins (or, equivalently, at least one of them) have 
zero X components and reach a contradiction. The initial as- 
sumption is used in an essential way to argue existence of the 
Bloch wall interpolating between spins with y components 
equal to +1 and —1. It is this part of the argument that we 
think is missing in Ref. |33l Note, that this argument applies 
to strong, as well as to weak random fields, so that the ground 
state is never, strictly speaking, field-dominated and always 
exhibits magnetization in the a;-direction. Moreover, the ar- 
gument does not depend on the dimension of the system, ap- 
plying in particular in one dimension. We argue below that in 
dimensions greater than one the effect still holds at small pos- 
itive temperatures, the critical temperature depending on the 
strength of the random field (and presumably going to zero as 
the strength of the field increases). 



D. RFIO at low positive T 

To study the system at low positive T, we need to ask what 
are the typical low energy excitations from the ground state. 
For e = 0, continuous symmetry allows Bloch walls, i.e. con- 
figurations in which the spins rotate gradually over a large 
region, for instance from left to right. The total excitation 
energy of a Bloch wall in 2D is of order one, and it is the pres- 
ence of such walls that underlies the absence of continuous 
symmetry breaking. However, for e > 0, a Bloch wall carries 
additional energy, coming from the change of the direction of 
the y component of the spin, which is proportional to the area 
of the wall (which is of the order L^ for a wall of linear size L 
in two dimensions), since the ground state spins are adapted 
to the field configuration, and hence overturning them will in- 
crease the energy per site. Similarly, vortex configurations, 
which are important low-energy excitations in the nonrandom 
XY model, are no longer energetically favored in the presence 
of a uniaxial random field. 

We are thus left, as possible excitations, with sharp domain 
walls, where the x component of the spin changes sign rapidly. 
To first approximation we consider excited configurations, in 
which spins take either their ground state values, or the reflec- 
tions of these values in the j/-axis. As in the standard Peierls 
argument^i, in the presence of the right boundary conditions, 
such configurations can be described in terms of contours 7 
(domain walls), separating spins with positive and negative x 
components. If m^ is the value of the x component of the spin 
ai in the ground state with the right boundary conditions, the 
energy of a domain wall is the sum of rairrij over the bonds 
[ij) crossing the boundary of the contour. Since changing the 
signs of the x components of the spins does not change the 
magnetic field contribution to the energy, the Peierls estimate 
shows that the probability of such a contour is bounded above 
by exp(-2/3X](ij) ™i™j): with (3 = J/ksT. 



We want to show that for a typical reahzation of the field, h, 
{i.e. with probability one), these probabilities are summable, 
i.e. their sum over all contours containing the origin in their 
interior is finite. It then follows that at a still lower T, this 
sum is small, and the Peierls argument proves that the system 
magnetizes (in fact, a simple additional argument shows that 
summability of the contour probabilities already implies the 
existence of spontaneous m). To show that a series of ran- 
dom variables is summable with probability one, it suffices to 
prove the summability of the series of the expected values. We 
present two arguments for the last statement to hold. 

If the random variables nii are bounded away from zero, 
i.e. rrii > ^/c, for some c > 0, the moment generating func- 
tion of the random variable J2iii) iTT'i'mj satisfies 



E 



expl 






< exp[-c/3L(7)], (3) 



with ^(7) denoting the length of the contour 7. The sum 
of the probabilities of the contours enclosing the origin is 
thus bounded by ^ exp[— c/3L(7)]. The standard Peierls- 
Griffiths bound proves the desired summability. 

The above argument does not apply if the distribu- 
tion of the ground state, m, contains zero in its sup- 
port. For unbounded distribution of the random field 
this may very well be the case, and then another argu- 
ment is needed. If we assume that the terms in the sum 



'YliUj) 'rrii'mj are independent and identically distributed, then 

E[exp(-2/3X;(y)W^mj)] = E[exp(-2/3m,mj)]^('^) = 
exp{L(7) logE[exp(— 2/3mimj)]} and we just need to ob- 



serve that E[exp(— 2/3mimj)] 



as /3 



cx) (since 



the expression under the expectation sign goes pointwise 
to zero and lies between and 1) to conclude that 
E[exp(— 2/3^,- J mimj)] behaves as eyi\)[—g{[3)L{-^)] for a 
positive function g{(3) with g{l3) ^ cx3 as /3 ^ 00. While 
mimj are not, strictly speaking, independent, it is natural to 
assume that their dependence is weak, i.e. their correlation 
decays fast with the distance of the corresponding bonds (ij). 
The behavior of the moment generating function of their sum 
is then qualitatively the same, with a renormalized rate func- 
tion g{(3), still diverging as /3 — > 00. As before, this is enough 
to carry out the Peierls-Griffiths estimate which implies spon- 
taneous magnetization in the a;-direction. We remark that our 
assumption about the fast decay of correlations implies that 
the sums of niimj over subsets of Z^ satisfy a large deviation 
principle analogous to that for sums of independent random 
variables and the above argument can be restated using this 
fact. 



E. Numerical Monte Carlo simulations 




Figure 2: (color online) Results of the Monte-Carlo simulation for 
the classical 2D-XY model in a 200 x 200 lattice. The Inset is a mag- 
netification of the main figure close to T = 0.7 J/ ka- 



power law decay of correlations, are particularly strong in the 
XY model in 2D. In particular, the 2D-XY model shows fi- 
nite magnetization (m) in small systems-'*, so that RFIO is 
expected to result in an increase of the magnetization. 

We have performed numerical Monte-Carlo simulations^ 
for the 2D-XY classical model [Hamiltonian ([T]i, with e = 1]. 
We generate a random magnetic field, h^ — rjiBy in the y di- 
rection. The r^i's are independent random real numbers, uni- 
formly distributed in [-x/SA/i^', VSA/i^']. Note that A/i'^ is 
thus the standard deviation of the random field h^. Bound- 
ary conditions on the outer square correspond to ctj = (1,0) 
[see Fig. [T]. The calculations were performed in 2D lattices 
with up to 200 X 200 lattice sites for various temperatures. The 
results are presented in Fig.|2] 

At very small temperature, the system magnetizes in the ab- 
sence of disorder (m approaches 1 when T tends to 0) due to 
the finite size of the lattice^*. In this regime, a random field 
in the y direction tends to induce a small local magnetization, 
parallel to h^, so that the magnetization in the x direction, m,, 
is slightly reduced. At higher temperatures (T ~ 0.7 J/ ks in 
Fig.lU, the magnetization is significantly smaller than 1 in the 
absence of disorder. This is due to non-negligible spin wave 
excitations. In the presence of small disorder, these excita- 
tions are suppressed due to the RFIO effect discussed in this 
paper. We indeed find that, at T = 0.7 J/k^, m increases by 
1 .6% in presence of the uniaxial disordered magnetic field. At 
larger temperatures, excitations, such as Bloch walls or vor- 
tices are important and no increase of the magnetization is 
found when applying a small random field in the y direction. 



III. RFIO IN OTHER SYSTEMS 

The RFIO effect predicted above may be generalized to 
other spin models, in particular those that have finite correla- 
tion length. Here we list the most spectacular generaUzations: 



Based on the above discussion it is expected that the RFIO 
effect predicted here will lead to the appearance of magneti- 
zation, TO, in the x direction of order 1 at low temperatures 
in systems much larger than the correlation length of typi- 
cal excitations. For small systems, however, the effect may 
be obscured by finite size effects, which, due to long-range 



A. 2D Heisenberg ferromagnet (HE) in random fields of 
various symmetries 

Here the interaction has the same form as in the XY case, 
but spins take values on a unit sphere. As for the XY Hamil- 



tonian, if the random field distribution has the same symmetry 
as the interaction part, i.e. if it is symmetric under rotations in 
three dimensions, the model has no spontaneous magnetiza- 
tion up to 4D (see Ref. l24ll25b . If the random field is uniaxial, 
e.g. oriented along the z axis, the system still has a contin- 
uous symmetry (rotations in the xy plane), and thus cannot 
have spontaneous magnetization in this plane. It cannot mag- 
netize in the z direction either, by the results of Ref. |25l Cu- 
riously enough, a field distribution with an intermediate sym- 
metry may lead to symmetry breaking. Namely, arguments 
fully analogous to the previous ones imply that if the random 
field takes values in the yz plane with a distribution invariant 
under rotations, the system will magnetize in the x direction. 
We are thus faced with the possibility that a planar field dis- 
tribution breaks the symmetry, while this is broken neither by 
a field with a spherically symmetric distribution nor by a uni- 
axial one. 



The Mermin-Wagner theorem.^^ implies that the model has no 
spontaneous magnetization at positive temperatures in 2D. For 
D > 2 spin wave analysis^-^dSiil shows the existence of spon- 
taneous magnetization (though a rigorous mathematical proof 
of this fact is still lacking). In general, one does not expect 
major differences between the behaviors of the two models at 
T y^ 0. It thus seems plausible that the presence of a ran- 
dom field in the quantum case is going to have effects simi- 
lar to those in the classical Heisenberg model. Similarly, one 
can consider the quantum Heisenberg antiferromagnet (HAF) 
and expect phenomena analogous to the classical case, despite 
the fact that unlike their classical counterparts, the quantum 
HF and HAF systems are no longer equivalent. We expect 
to observe spontaneous staggered magnetization in a random 
uniaxial XY model, or random planar field HF. A possibility 
that a random field in the z-direction can enhance the anti- 
ferromagnetic order in the xy plane has been pointed out in 
Ref.Hl 



B. 3D XF and Heisenberg models in a random field of various 
symmetries 

We have argued that the 2D XY model with a small uni- 
axial random field orders at low T. Since in the absence of 
the random field spontaneous magnetization occurs only at 
T = 0, this can be equivalently stated by saying that a small 
uniaxial random field raises the critical temperature Tc of the 
system. By analogy, one can expect that the (nonzero) Tc of 
the XY model in 3D becomes higher and comparable to that 
of the 3D Ising model, in the presence of a small uniaxial field. 
A simple meanfield estimate suggests that Tc might increase 
by a factor of 2. The analogous estimates for the Heisenberg 
model in 3D suggest an increase of Tc by a factor 3/2 (or 3) 
in a small uniaxial (or planar rotationally symmetric) field re- 
spectively. These conjectures are the subject of a forthcoming 
project. 



C. Antiferromagnetic systems 

By flipping every second spin, the classical ferromagnetic 
models are equivalent to antiferromagnetic ones (on bipartite 
lattices). This equivalence persists in the presence of a ran- 
dom field with a distribution symmetric with respect to the 
origin. Thus the above discussion of the impact of random 
fields on continuous symmetry breaking in classical ferromag- 
netic models translates case by case to the antiferromagnetic 
case. 



IV. TOWARDS THE EXPERIMENTAL REALIZATION OF 
RFIO IN ULTRACOLD ATOMIC SYSTEMS 

Further understanding of the phenomena described in this 
paper will benefit from experimental realizations and inves- 
tigations of the above-mentioned models. Below, we dis- 
cuss the possibilities to design quantum simulators for these 
quantum spin systems using ultracold atoms in optical lattices 
(OL). 



A. Two-component lattice Bose gas 

Consider a two-component Bose gas confined in an OL 
with on-site inhomogeneities. The two components corre- 
spond here to two internal states of the same atom. The 
low-T physics is captured by the Bose-Bose Hubbard model 
(BBH)"*-' (for a review of ultracold lattice gases see Ref. 1441): 



Hbbh - E [y"^K - 1) + ^N,(N, - 1) 
j 



(4) 



5][(jbb]b, + JBB]B,)+h.c. 



n 



-E ?btB,+h.c, 



D. Quantum systems 

All of the effects predicted above should, in principle have 
quantum analogs. Quantum fluctuations might, however, de- 
stroy the long-range order, so each of the discussed mod- 
els should be carefully reconsidered in the quantum case. 
Some models, such as the quantum spin S — 1/2 Heisenberg 
model, for instance, have been widely studied in literature^^. 



where hj and B^ are the annihilation operators for both types 
of Bosons in the lattice site j, rij = bjbj and Nj = BJBj are 
the corresponding number operators, and {j,l) denote a pair of 
adjacent sites in the OL. In Hamiltonian (01, (i) the first term 
describes on-site interactions, including interaction between 
Bosons of different types, with energies Uh, C/r and U^b', (ii) 
the second accounts for on-site energies; (iii) the third de- 
scribes quantum tunneling between adjacent sites and (iv) the 




Figure 3: (color online) Atomic level scheme of a two-component 
Bose mixture in a random optical lattice used to design spin models 
in random magnetic fields (see text). 
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fourth transforms one Boson type into the other with a prob- 
ability amplitude |il|/?l. The last term can be implemented 
with an optical two-photon Raman process if the two Bosonic 
'species' correspond to two internal states of the same atom 
(see also Fig. [3]). Possibly, both on-site energies Vj, Vj and 
the Raman complex amplitude ilj can be made site-dependent 
using speckle laser light' ^'^"^-'^ . 



Consider the limit of strong repulsive interactions (0 < 
Jb, JbjI^jI ^ Uh,UB,UBb) and a total filling factor of 1 
(i.e. the total number of particles equals the number of lat- 
tice sites). Proceeding as in the case of Fermi-Bose mixtures, 
recently analyzed by two of the authors in Ref. S, we derive 
an effective Hamiltonian, i/ett, for the Bose-Bose mixture. In 
brief, we restrict the Hilbert space to a subspace Eq generated 
by { Jl \nj ,Nj)} with rij+Nj = 1 at each lattice site, and we 
incorporate the tunneling terms via second-order perturbation 
theory as in Ref. |23i We then end up with 



Y 



E^.-^.-E V^.+h-c 



(5) 



where S, = W-HjV, V is the projector onto Eq and Mj = 



B -Bj. Hamiltonian i/etf contains (i) a hopping term, Jji, (ii) 
an interaction term between neighbour sites, Kj,i, (iii) inho- 
mogeneities, Vj, and (iv) a creation/annihilation term. Note 
that the total number of composites is not conserved except for 
a vanishing il. The coupling parameters in Hamiltonian (|5]l 



where dj,i = Vj — vi and A^ ; — Vj — Vi. Hamiltonian H^ff 
describes the dynamics of composite particles whose annihi- 
lation operator at site j is Bj — WaBjV. In contrast to the case 
of Fermi-Bose mixtures discussed in Ref. |23l where the com- 
posites are fermions, in the present case of Bose-Bose mix- 
tures, they are composite Schwinger Bosons made of one B 
boson and one h hole. 

Since the commutation relations of Bj and B] are those 
of Schwinger Bosons"^^, we can directly turn to the spin 
representation^ by defining SJ + iS^ = Bj and S| — 

1/2 — Mj, where Mj = B^Bj. It is important to note 
that since Raman processes can convert h Bosons into B 
Bosons (and conversely), J^ji-^j) '■^ not fixed by the total 
number of Bosons of each species, i.e. the z component of 
TO, J^ji^j) is 11°'^ constrained. For small inhomogeneities 
i^jd = ^3 ~ ""h^jd = V, - V; < Ub, Ub, Ubs), Hamiltonian 
i/eff is then equivalent to the anisotropic Heisenberg XX Z 
model"*^ in a random field: 



H, 



eff 





Sf + S] 


sf)- 
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E ^^ 
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+ /.|S| 
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3) 



where 



J± 



4JbJB 



C/bB 



(9) 



(10) 



J. = 2 



2Jl 

Ub 



2_Jl 
Ub 



■^b 
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"-3 



n^. 



^ 



-n\ 



UbB 

; /i? 



V, 



(11) 
CJz/2 , (12) 



' The coupling parameters are the same as calculated in Refs. 1 13231 except 
for the third term in Eq. I?) which corresponds to a virtual state with two 
B bosons in the same lattice site — forbidden for Fermions. 



with (^ the lattice coordination number, V, 



V, 



C[4j2M + 4j2/C/B-(j2 + J|)/f/bB] and% - n^^ + i^]. In 
atomic systems, all these (possibly site-dependent) terms can 
be controlled almost at will^Mii^. In particular, by employing 
various possible control tools one can reach the HF {J± = J^) 
and XY {Jz = 0) cases. 



B. Bose lattice gas embedded In a BEC 

The quantum ferromagnetic XY model in a random field 
may be alternatively obtained using the same BBH model, but 
with strong state dependence of the optical dipole forces. One 
can imagine a situation in which one-component (say b) is in 
the strong interaction limit, so that only one b atom at a site is 
possible, whereas the other (B) component is in the Bose con- 
densed state and provides only a coherent BEC 'background' 
for the fe-atoms. Mathematically speaking, this situation is de- 
scribed by Eq. (|4|i, in which n^'s can be equal to or 1 only, 
whereas Bi's can be replaced by a classical complex field 
(condensate wave function). In this limit the spin S — 1/2 
states can be associated with the presence, or absence of a b- 
atom in a given site. In this way, setting Vj = and fl' = 0, 
one obtains the quantum version of the XY model ([Tji with 
J = Jb and a uniaxial random field in the x direction with the 
strength determined by fl^. 



C. Two-component Fermi lattice gas 

Finally, the S — 1/2 antiferromagnetic Heisenberg model 
may be realized with a Fermi-Fermi mixture at half filling 
for each component. This implementation might be of spe- 
cial importance for future experiments with Lithium atoms. 
As recently calculated"*^, the critical temperature for the Neel 
state in a 3D cubic lattice is of the order of 30nK. It is well 
known that in a 3D cubic lattice the critical temperatures for 
the antiferromagnetic Heisenberg, the XY and the Ising mod- 



els are Tc.xy ^ l.bTc.Heis, and Tcjs 



2Tr 



He 



The 



estimates of these critical temperatures can be, for instance, 
obtained applying the Curie-Weiss mean field method to the 
classical models. Suppose that we put the Heisenberg anti- 
ferromagnet in a uniaxial (respectively, planar) random field, 
created using the same methods as discussed above, i.e. we 
break the SU{2) symmetry and put the system into the uni- 
versality class of Xy (respectively, Ising) models. Mean field 
estimates suggest then that we should expect the increase of 
the critical temperature by factor 1 .5 (respectively, 2), that is 
up to ~ 45 (respectively, 90)nK. Even if these estimates are 
too optimistic, and the effect is two, three times smaller, one 
should stress, that even an increase by, say lOnK, is of great 
experimental relevance and could be decisive for achieving of 
antiferromagnetic state. 

We would like to stress that similar proposals, as the three 
discussed above, have been formulated before ''^, but none of 
them treat simultaneously essential aspects for the present 
schemes: i) disordered fields, but not bonds; ii) arbitrary di- 
rections of the fields; iii) possibility of exploring Ising, XY 



or Heisenberg symmetries; iv) realizing the coherent source 
of atoms; and v) avoiding constraints on the magnetization 
along the z axis. 

It is also worth commenting on what are the most impor- 
tant experimental challenges that have to be addressed in or- 
der to achieve RFIO. Evidently, for the proposals involving 
the strong interaction limit of two-component Bose, or Fermi 
systems, the main issue is the temperature which has to be of 
order of tens of nano-Kelvins. Such temperatures are start- 
ing to be achievable nowadays (for a carefull discussion in the 
context of Fermi-Bose mixtures see Ref. 49), and there exist 
several proposals for supplementary cooling of lattice gases, 
using laser (photons) or couplings to ultracold BEC (phonon 
cooling) that can help (for reviews see Ref. l44M6l 



V. SUMMARY 

In this paper, we have proposed a general mechanism 
of random-field-induced order (RFIO), occurring in systems 
with continuous symmetry, placed in a random field that 
breaks, or reduces this symmetry. We have presented rigor- 
ous results for the case of the 2D-classical ferromagnetic XY 
model in a random uniaxial field, and proved that the system 
has spontaneous magnetization at temperature T — 0. We 
have presented also a rather strong evidence that this is also 
the case for small T > 0. Several generalizations of this 
mechanism to various classical and quantum systems were 
discussed. We have presented also detailed proposals to re- 
alize RFIO in experiments using two-component Bose lattice 
gases, one-component Bose lattice gases embedded in BECs, 
or two-component Fermi lattice gases. Our results shed light 
on controversies in existing literature, and open the way to 
realize RFIO with ultracold atoms in an optical lattice. 

It is worth mentioning two further realizations of RFIO 
studied by us recently. RFIO occurs in a two-component 
trapped Bose gas at T = 0, when the gas is condensed and 
the two components are coupled by Raman transition of ran- 
dom strength, but fixed phase. Although such a system be- 
longs to the universality class of the (trapped, i.e. located in 
an inhomogenous field) XY model, it exhibits the RFIO ef- 
fect in a much stronger manner than the XY model discussed 
in the present paper We have found this observation important 
enough to devote a separate detailed paper to it^*'. Similarly, 
we have studied numerically RFIO in ID for quantum XY 
and Heisenberg chaina^i. In such systems, even at T = 0, 
magnetization vanishes, but amazingly enough the RFIO ef- 
fect seems to work at the level of the magnetic susceptibilities. 
Adding a random field confined to a certain axis (respectively, 
plane), increases siginificantly the magnetic suceptibility in 
the perpedicular directions (respectively, direction). 
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